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About distribution natural waves in dissipative 
inhomogeneous flat bodies 
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In this article propagation of natural waves in dissipative inhomogeneous planar bodies is considered. Wave motions are described by linear 
integral-differential equations. Solving this problem, we obtain a relationship between the wave velocity and its length. The task of this kind is 
of great interest for geophysicists, in the field of engineering and construction. 


INTRODUCTION 

Many constructions and engineering structures, working in dynamic conditions, consist deformable bodies with different viscoelastic properties 
[1,2,3,4,6,7]. In addition, wave processes in elastic bodies play an important role in connection with the processing of signals, in particular in 
connection with the creation of mechanical resonators and filters [8, 9, 10]. The mechanisms by which the energy of elastic waves is converted into 
heat are not entirely clear. Various loss mechanisms are proposed [11 + 16], but not one of them does not fully meet all the requirements. Probably 
the most important mechanisms are internal friction in the form of sliding friction (or sticking, and then slipping) and viscous losses in pore fluids; 
the latter mechanism is most significant in strongly permeable rocks. Other effects that are probably generally less significant are the loss of some of 
the heat generated in the phase of compression of wave motion by thermal conductivity, piezoelectric and thermoelectric effects and the energy going 
to the formation of new surfaces (which plays an important role only near the source). Therefore, the development of a unified methodology and 
calculation algorithm, wave fields of dissipative inhomogeneous layered bodies, is a topical problem of the mechanics of a deformable solid [17, 18, 
19]. 


Formulation of the problem 
Let the Cartesian (x, y, z) coordinate system, with the origin and the OZ axis, be given a sequence of parallel planes (Fig. 1) 


y = hithoet...+hn 


y = hithot...+hn-1 


Ant 2 Mv 2 Poa 


Figure 1 The design scheme: the body on the half-space 
Z=—O,2 Sh, 2h hy + ex. Lo = hi, hy Po Hi, 


Plane Z = h, + i; Paarl , (N =2) will be called 1—Mhorizon Suppose that the spaces between the planes mentioned are filled with 
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isotropic elastic media forming parallel layers. Layers 0<Z <h, characterized by permanent Ais Ho; Po, , will be called zero. Wednesday, 


however, h, =F h, Hea h,, <Z< h, =F He Esa ce i, +h, filling the space between —Mandn + 1—mbhorizons characterized by 


n+l 


constant Ay »H,,P 1, + will be called 2 — M1 layer. It will always be assumed that adjacent layers differ from each other in at least one of the 


constants A e) Lt andp. In the theoretical study of the described processes, we shall assume that within each layer the wave propagation is described 


by the usual equations of the theory of elasticity. As for the conditions on the interfaces of adjacent layers, we assume that the components of the 
vector of elastic displacements and the stress tensor remain continuous when passing through them. This contact is called hard [20]. The dynamics of 
dissipative inhomogeneous two-layer flat structures is investigated in the article. 

Accounting for internal friction, caused by the dissipation of energy in the material of structures, is a more difficult task. Soft layers of multilayer 
structures (aggregates), as a rule, are made of materials with developed rheological properties. Therefore, the scattering of energy must first of all be 
taken into account for soft layers, since it mainly occurs precisely when deforming these layers. Mechanical systems, for which the viscoelastic 
properties of their elements are identical, will be called dissipative homogeneous, systems with different rheological characteristics — dissipative 
heterogeneous [21,22,23]. 


Equations of motion of the deformed layer in the absence of mass forces have the form [9]: 


2- 
~ 2- a ~ . > O Uu : 
HVut(A,+ HL, ygrad di9 u= P; Ar? ,(G=1,2,3..) () 
Here u(u yo y? u ‘) -vector of displacement of points of the medium; - the density of the material; ui — displacement components; v Jj - 


Poisson's ratio; 


7 VE; ~ VE; 


Ani 3 sa 
/~ d+v,)0—2v,) cd 2d+v,)° 


~ 


E- Operational modulus of elasticity, which have the for [9,13]: 


E,o()= E,) olt)—| Ry (trol |v 


0 


gt )- arbitrary time function; R Ej (t =— )-the core of relaxation; Bx —instantaneous modulus of elasticity; We assume the integral terms in (2) to 


be small, then the functions gt ) =VWV (t Je ‘eet , where YW (t )-a slowly varying function of time, @p -real constant. Further, applying the freezing 


procedure [9], we note relations (2) with approximations of the form 


E,p=E,,|I-T¢ (,)-i} (@)p. 


where “(o, )=[R, (r)cos @,t drt. I? (@,)= 


0 


R, (c)sin Opt dT. respectively, the cosine and sine Fourier images of the relaxation 


ota 8 


_ “Bit 741-4; 
core of the material. As an example of a viscoelastic material, we take three parametric relaxation nuclei R j t)=A ie /t . On the 


influence function R j (t —_ t) The usual requirements of inerrability, continuity (besidesé = T), signs - certainty and monotony: 
Dr 
dR(t) 
R > 0, —— <0, 0< ] R(tHdt<1. 
dt 
0 


Ul - vector of displacements of the medium j-the layer. 
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Two types of conditions can be defined on the boundary of two bodies: 


1. In the case of a rigid contact at the interface, the condition of continuity of the corresponding components of the stress tensor and displacement 
vector is set, i.e. 


Oo. 2. -.0 2.0), 
Oy, =O, 3 On =On 3 


yy xy xy 
aye Oa) (3a) 
uy =U, ; uy =uy’. 


If there is no friction at the interface, 


() _ (2), () _ (2) _¢. ( _,,(2) 
Gy Oy Oy ~Oy HUs Wy H's (2,6) 
dd) _p. (1) _ 
2. On the free surface, the condition of freedom from stress is set,r.e. O = 0 ’ C |. ~ 0, 
(2,c) 


whave 
Ou, . ou, Og. 
On =4,0,+2u, oy) =n) ——+— 


Oy ox 


D729 42 j Ou; 08; 
eal As ae ae 


ox oy 


Methods of solution 


Now consider the solution of the differential equation (1) - (2) for one layer. The equation of motion in displacements reduces to the following form: 


_( 07u, u'r + 4+@) 0(Ou, ; 03,) Ou, =i 
rae aye | ele Oy) ee 
(3) 

_(O9 Ar) ~ _.dfdu, dS, Og, 

Ela 5 (+A, +4.) 2 Pio 

Ox Oy Oy\ Ox Oy Ot 

h P, - density of the material. We find the solution of the problem in the form: 
where 

m=aU0 er’ 8 ave? n=l2aaaN (4) 


h U ny) ) and V, (y) — amplitude complex vector - function; k is the are number; C = C Rt iC ; ~ complex phase velocity; ao — complex 
where 


frequency. 
To clarify their physical meaning, consider two cases: 
1) k = kr; C = Cr +iC,i, then the solution (4) has the form of a sinusoid with respect to x, the amplitude of which decays in time; 


2) k= kr +iki; C = Cr, Then at each point x the oscillations are steady, but with respect to x they decay. 
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In both cases, the imaginary partskiorCicharacterize the intensity of dissipative processes. Substituting (4) no (3), we obtain: 


— 2 ar = ae Ae 2772 
Z,(U" -KU, )+ (4, +H, ) ik(KU, +V)+ p,17C'U, =05 65, 
— 2 a — \.7(/. 2 72 
7,(U" -kU, )+(A, +2, ) ik(ikU, +V!)+ p,PCU, =0. 
Thus, we have equations (5) of the second order for two domains each. We solve the problem directly, without reducing the equation to a fourth- 
order equation. All the arguments are given for the layer. 


We find the particular solution of the system (5) in the form 


U A 


n n “ry 


V B 


n n 


‘i tm — contact. A homogeneous algebraic system with respect toAnandBuhas non-trivial solutions if its determinant is zero 
where 


CPSC 20 aCe (C2, -C2, ker, 
= 0, (6) 
i(Cy, —Ci,, kr, C0 a a1C7, Sea +C°k? 


c = (A, + 2H)! P,- es = H,, 1 Des atn = 0 magnitudes 


are respectively the velocity of the waves of compression and shear in an 
elastic medium [4]. Equation (6) can have four roots 


Care =H 1=-C"/ Ce: (7,)3.4 St —C iC, n—01. a 4. 


As a result, we find four particular solutions of the form 


4 

U,, = C 4, (7 i” 

y => nil p e »n=0,1. (7) 
n i=1 


AL 


Substituting the values (rn) iB (7) we findAni, Bniattn =(1n) i. 


Expressions for the offset are as follows: 


U, ik ae ik e& 
=C,, . | en Ch oki 4 
Vv, — kq, kq, 
C kS, —ksyy C kS, +ksyy , 
+™13 14 > 
—ik —ik 
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Consequently, for both the hard and sliding contacts, we obtain a set of six boundary conditions that lead to six homogeneous equations with six 
unknownsC 11 ,C12,C13 ,C14 ,C22, C24. For such a system of equations to have non-trivial solutions, the determinant of the coefficients must be zero. The 
last equation gives the dispersion equation for dissipative systems, 


where 5, =U-C*7/C7n)'?; G, =A-C7/C’n.), n=0,1 


As an example, let us consider the problem of propagation of natural waves in a viscoelastic layer on a half-space. 


Hard contact. The dispersion equation has the following form 


(1+ 57) een (1+ 57) eet _ 2S 28% O O 
_s _ se _ 1 —eq _ 1 
— 29,7" 24a," [s + +}. [s + + ei O O 
Sy Sy 
(i+s7)e% .(1+57)e% 2% aan (Ls 577, -—27y, 
=0 
; 1)2z 1 2g, 
a 2g; ody 2G," [s + + Jee [s zs Le "q1 f \ty 
Sy Sy v1 > 
em em ; em em —] a | 
= _- & les les _ 1 
gqe™ 7 se ae” seg? —q = 
s Ss S 
(8) 


where ¢ — dimensionless wave number = kh,yi= Ly, / o orA+2 = 


Ae’ 
= 2 (1-v)/(1-2v). . As the relaxation nucleus of a viscoelastic material, we take a three-parameter core R(t ) = 


Rizhanitena-Koltunova [13], 


which has a weak singularity, where A, @, B -parameters materials [13]. We take the following parameters: A = 0,048; Bp =0,05; a=0,l 


Using the complex representation for the elastic modulus, described earlier. 


The roots of the frequency equation are solved by the Mueller method, at each iteration of the Muller method is applied by the Gauss method 
with the separation of the main element. Thus, the solution of equation (8) does not require the disclosure of the determinant. As the initial 


approximation, we choose the phase velocities of the waves of the elastic system. For free waves with R j =0 the phase velocities and the wave 


number are real quantities. No the calculations, we acce the following parameters values: 
0 = pi/p2 =0, 75; B = 104 sn = 1. 


Let us consider two variants of the dissipative system. In the first variant, the dissipative system is structurally homogeneous. 

Wave number & varies from 0 to 3. The results of the calculations are shown in Fig. 2a. Dependence of frequencies and damping on the 
dimensionless wave numberéturned out to be monotonous, and the character of the dependence is the same for the frequencies and damping 
coefficients. In the second variant, the dissipative system is structurally inhomogeneous: the half-space under consideration, equation (8), and elastic 
parameters coincide with those adopted above. The results of the calculations are shown in Fig. Dependence of frequencies onEThe same as for a 
homogeneous system: the corresponding curves coincide with an accuracy of up to 5%. Dependence of damping factors on€no monotonic. 
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Of particular interest is the minimum valueéwith a fixed damping factor: 


d=min (-ax ), k= 1,2,.....K, (9) 


hered — coefficient that determines the damping properties of the system (we call it the global damping factor). 


Figure 2a Variation of complex Eigen frequencies from the wave number. A) Dissipatively homogeneous mechanical system 


For a homogeneous system, the coefficientdis entirely determined by the imaginary part of the first complex-frequency modulus. For an 
inhomogeneous system in the role of the coefficientéthe imaginary parts of both the first and second frequencies may act as a function of their values. 
"Change of roles" occurs with a characteristic valueéat this value, the real parts of the first and second frequencies are the closest. Coefficientdat the 
indicated characteristic value has a pronounced maximum. 


Sliding contact. The dispersion equation is similar in form to equation (8). All parameter values are the same as those used above. Figures 3a and 6b 


show the dependences of the frequencies and damping coefficients on the wave numberé, respectively for structurally homogeneous and 
inhomogeneous systems. The obtained results confirm the earlier conclusions. Change the parameter, from 
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1,0 


0,1 


Figure 2b Change in complex Eigen frequencies from the wave number. b) Dissipative inhomogeneous mechanical system 


which depends so significantly on the coefficient 5, can be achieved by varying the geometric dimensions of the elements without changing their 
mechanical properties. 


Figure 3b Dissipative homogeneous mechanical system 
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Figure 3b Dissipative inhomogeneous system 
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Figure 4 Dispersion curves of phase velocities KR ; XR? and attenuation coefficient An ; X12 for the duralumin strip in contact with the 


acrylic medium 
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Thus, a promising opportunity for effective control of the damping characteristics of structurally inhomogeneous viscoelastic systems opens up 


by changing their inhomogeneous systems with close frequencies. As a second example, let us consider the propagation of natural waves in a plane 
layer located in a deformable (viscoelastic) medium (Fig. 4). 

The results of the calculation are shown in Fig. 4. Dependence of frequencies on turned out to be the same as for a dissipative homogeneous 
system: the corresponding curves coincide with an accuracy of up to 5%. As for the coefficients of damping, their behavior has changed radically: the 
dependence wi~ € became no monotonic. The global damping factor at the specified characteristic value€ has a pronounced maximum. 


CONCLUSION 

1) In the course of solving the problem, the propagation of waves in dissipative - inhomogeneous media revealed no monotonic dependences of the 
damping rate on the physic - mechanical and geometric parameters of the system. In Dissipatively inhomogeneous media, the dependence of the 
phase velocity and the damping rate on the geometric and physic-mechanical parameters of the system turned out to be no monotonic; 


Cp1=5400 m/e; V 1=0.35;Cp2=2300 m/c; V =0.35; Cs1=1311 m/c; p1=0.138 xr/m?; Cy=1311 m/c; p2=0.126xr/M?. 


2) Based on the obtained numerical results, it is revealed that the possibility of detachment of thin-walled structures from a soft layer and the effect of 
magnitude up to the resonance speed on the dimensions of the contact area. Also taking into account the viscous properties of the material, 15 - 10% 
increases the values of the phase velocities; 

3) Revealed that the phase higher forms of the expansion and torsion waves exceed the highest possible speed (C) of waves in an infinite medium, the 
group velocity never exceeds C. Also found that the group velocity of 10-15% exceeds the non-dispersive medium, comparison by a dispersive 
medium. In other words, the forms of the pulses do not remain unchanged, as in homogeneous elastic bodies. 
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